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Abstract 

For the low energy Standard Model neutrino-matter interactions, we calculate neutrino 
pair (vv) emissivites in superfluid quark matter. Just below the critical temperature, 
Cooper pairs continuously break and recombine, resulting in the emission of vv pairs 
with a rate that greatly exceeds the standard quark modified Urea and bremsstrahlung 
rates. At the same baryon density in baryonic and quark matter, the ratio of baryon to 
quark vv emissivites lies in the range 2 — 5 for the densities of interest in the long-term 
cooling of solar mass compact stars. We also find that in matter containing hyperons, 
vv emission can occur with hyperons of all species. 

PACS: 97.60.Jd, 24.85.+p, 26.60,+c, 95.30.Cq 

In a prescient paper, Flowers, Ruderman, and Sutherland |[| showed that just below the critical 
temperature T c , the neutrino pair (vv) emissivity from the breaking and recombination of superfluid 



neutron pairs (n + n ^n + n + v + v) greatly exceeds that from the so-called modified Urea processes 
(n + n — >n + p + e~ + v and n + p — > n + n + e + + v) in neutron star interiors. It is only recently that 
the pair breaking and formation (PBF) processes have been included in calculations of the long term 
cooling of neutron stars M. The principal effect of the PBF process is to rapidly cool the star during its 
. 5^ early hundreds of years of thermal evolution, until competing //-emitting processes (for recent reviews 
see Ref. [Q]) begin to dominate the further cooling of the star for several thousands or million years that 
they might remain observable with such instruments as HST, Chandra, and XMM. The precise extent 
to which the core temperature falls depends on the smallest of the pairing gaps among the superfluid 
species present in the star (cf. and references therein). To date, the role of PBF processes has been 
explored only for baryon superfluidity @. 

In this work, we calculate the vv emissivity in superfluid quark matter likely to be encountered in 
neutron star interiors ||. For the low energy Standard Model neutrino- matter interactions, we show 
that the predominant contribution to the total vv emissivity comes from the vector channel, while the 
axial channel provides a small correction. This heirarchy is similar to that found in baryonic matter [|l| . 
We also consider the effect of quark masses on the emissivities in limiting cases. We compare our results 
for vv emissivities with those of the modified Urea and bremsstrahlung processes in both quark and 
nucleonic matter @-f|. Expected vv emissivities from hyperons in baryonic matter and from p-wave 
pairing in quark matter are also addressed. 
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The Hamiltonian that describes the neutral current quark-neutrino weak interaction is given by || 



G 

H W = ^= W'^icv - C A j5)ll>i,a] [^7a(1 - 7s)V^] , (1) 

here Gp is the Fermi constant for the weak interaction, and cy and ca are flavor dependent vector 
and axial-vector coupling constants, respectively. The index i denotes quark flavor (Nf = 3) and a 
denotes color (iV c = 3). The index i stands for lepton flavor e, fi or r. At leading order (tree level), 
flavor changing processes do not occur through the weak interaction, hence the diagonal structure in 
flavor and color space. Unlike baryons, quarks in a superfluid phase can lock themselves in preferred 
directions in color and flavor space ||. However, we are considering neutrino pair {yv) production 
mediated by the neutral current weak interaction, which is blind to color and flavor. These degrees of 
freedom are, however, counted in the density of states at the Fermi surface. 

We begin with the case in which quarks of all flavors are taken to be massless. In this case, the full 
structure of the weak interaction through 4-component Weyl spinors is needed (In the nonrelativistic 
case, e.g. for baryons in the superfluid phase, the V-A interaction has large matrix elements only in 
the 2x2 block diagonals in the space of 4-component spinors. This allows for simpflication in that 
2-component spinors with a suitable normalization are adequate.) The helicity basis (which is also the 
chirality basis for massless particles when we identify positive helicity with right handed particles {R} 
and negative helicity with left handed {L} ones) is convenient in this case. The quark field operator 
may be expressed as [] 



= E ( e iTi -%{Xp,Rh,R + Xp,iA>,l} + e l »-\{Xp,nbl iL - Xp,Lbl R } J , (2) 
p \ ' 

where x's are 4-component Weyl spinors and we have paired particles of opposite helicity according to 

bp,R = u p a$ iR + ^p a -p,L and b P,L = u p a p,L ~ v p a -p, R ■ (3) 

The quasiparticle excitations are determined through a Bogoliubov-Valatin transformation on the 
usual creation and annihilation operators (a and a'), and are distributed according to the Fermi-Dirac 
function at finite temperature. For isotropic 1 So pairing, the coherence factors u v and v p satisfy 

A ^ = ifl + #\ ^d «2 = i(l-£V (4) 



1Ep p 2 V EpJ p 2 V Ep 



The quasiparticle excitations carry momentum p and energy E p = ye p 2 + A 2 , where e p = p—fi and A is 
the superfluid gap. The field operator in Eq. (0) acts on a quasiparticle state \p R , p' L ) and restores it to 
the condensate. Conversely, it can also excite the condensate so that pair breaking and recombination 
occur simultaneously. We denote the matrix element for this process by M q and the neutrino matrix 
element by M v . The time dependence on the field operator has been dropped, since the quasiparticle 



In this form, the quark field operator does not take into account the fact that diquarks may condense in 
preferred channels in color and flavor space, as, for instance in the three flavor color-flavor locked (CFL) phase, 
or the two flavor superconducting (2SC) phase (see, for example, [0]). We return later to comment on the 
expected modifications in "QCD condensed matter" . 
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lifetime is essentially infinite compared to the weak interaction time. Then, the expression for the 
emissivity is 



e - l N N n( Gf \ 2 ^ )A f d3pd3p ' f(F)f(F 
e--N v N } NA- M \ —zj ——-f{E p )f{E pl 



\V2J (2vr) 12 7 2 Po 2p' 

r pipl {qi0 + q2o) \ Mq \ 2 \M^ 6\p + p' - gi - q 2 ) , (5) 
J 2q w 2q 20 

where p and p' denote the four-momenta of quarks and q\ and q 2 those of the neutrinos. The number 
of neutrino flavors is denoted by N v . The simple factor of NfN c outside arises, because we have 
disregarded strong interaction effects and flavor changing processes that occur at higher orders. The 
overall factor of 1/2 avoids double counting of the same quasiparticle collision. / denotes the Fermi- 
Dirac distribution function. Summing over the neutrino spins yields the neutrino squared matrix 
element, \M U \ 2 = S(q 2 ^q\ u + qi^u - fiV<?i-<?2 - i^apQiaQ2p)- The quark matrix element receives 
contributions from the vector (V), axial (A) and the mixed vector-axial (M) channels. 
The vector part of the interaction acting on the state \pr,p l ) gives 

i\ M q\ 2 )v = c 2 v {u p v p , + ^v) 2 (Xp,L7%-/,LXp', L 7 iy Xp,L) , (6) 



where we have used the spinor identity Xp,Rl a Xp' r — Xpi £7 a Xp,L |TlJ in writing the term in the 
rightmost bracket. Using the standard normalization for massless spinors, the vector contribution to 
\M q \ 2 becomes 



(\M q \ 2 ) v v = 2(u p v p/ + v p u p/ ) 2 ■ (7) 

where T = Tr [^(1 - 7 5 )^V(1 - 75)]. Utilizing the relation Xp' tR l a l5Xp,R = -Xp,Ll a l5Xp> tL = 
Xp,Ll a Xp> li ^ ne sc i uare d matrix elements for the axial and mixed contributions are found to be 



(|M ff | 2 )JST = 2(«pty - «ptv) a ' 4^ and (\M q \ 2 )^ = 2(u 2 pl v p 2 -v 2 pl u p 2 )-^^T. (8) 



The total contribution from the quark squared matrix element is simply the sum 

\M q \ 2 = \M q \y + |M g |^ + |M ? ||^. The contraction of the neutrino and quark squared matrix elements 

yields, up to differences in couplings and coherence factors, a generic expression of the form 

\M q \ 2 c x \M V \ 2 ex [16( 9l -p){q 2 -p') +48(g 2 ■ p)( qi ■ p')} , (9) 

where C denotes any of the channels V, A, and M. 

We turn now to the phase space integrations. The neutrino phase space integral is of the form 

N a p = [ ^_Sl ( L^l ( qia q 2f3 + q 2a q lf3 )(q w + q 20 ) 5 4 (qi + q 2 - Q) 
J IqiQ Zq 20 

= ^(Q 2 g a ? + 2QaQpMQo)®(Ql - IQI 2 ) , (10) 

— * 

where an additional delta function in the total neutrino energy Qq and momentum Q has been intro- 
duced. The obvious symmetry between a and f3 in the above integral may be exploited to simplify the 
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calculation (e.g. the €p Va p terms from the neutrino squared matrix element vanish once the neutrino 
integral is performed). The momentum conserving delta function in Eq. (||) can be used to eliminate 
the integral over p'. Next, we shift the p integral using p — > p + ^, and in view of the restriction 
that |Q| < Qq ~ A and A <C |p| ~ /i, the energy conserving delta function becomes S(Q — 2E P ). In 
addition, the two Fermi-Dirac functions both become f(E p ). | 

In the vector channel, the momentum shift p — > p + y allows us to write 

(u p v Q - p + v p u Q ^) 2 P ^ 2 u 2 p+ cLvl p+ z + v 2 p+ zu 2 _ p+ z + 2u p+ 2 v p+ g. v_ p+ zu_ p+ g_ . (11) 

Using Eq. (|j) and the expression for the quasiparticle energy, the first two terms become 

1A 2 xA 2 CjL , t ^\ t 



where x = A//i <C 1, justifying the neglect of the 0(x 2 ) terms. The third term becomes A 2 /(2E 2 ). 
Thus, the coherence factors for the vector channel take the form 

A 2 

(u p Vp> + v p u p/ ) 2 = — . (13) 
E p 

For the axial contribution to the quark matrix element, we retain terms up to 0(x 2 ) to get 

2 p^p+# (Q.p) 2 A 2 
[UpVQ-p - VpU Q -p) — > = — — — . (14) 

Eqs. ([IB] ) and (III]) enable the coherence factor for the mixed contribution to be written as 



/ 2 2 2 2s 
{u p ,Vp - V p ,Up ) 



A 2 (Q.p) 2 A 2 Q.pA 2 



(15) 



The mixed contribution, being odd in Q, vanishes on angular integration. Integrals over Q can be 
performed using the formula 

J dn q d\Q\ |Qp (|Q.p|)« = -i^-|p|«| d|Q| IQI-+-. (16) 

The integrals over Q Q , Q, and dQ p are now straightforward, leaving only the integral over |p|. We can 
make a change of variable at the Fermi surface from |p| to E p to cast the final result as a 1-dimensional 
integration over E p . Collecting these results and summing the emissivities in the vector (V) and axial 
(A) channels, the vv emissivity from PBF in matter comprised of massless quarks is 



2 It is possible to calculate corrections resulting from E p / E p /, but lead to more complicated expressions 



than found here (e.g, see Ref. p^ , |i3| ) 
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F{y) 



1 fCA 



7 VCy 



(17) 



where z/(0) = NfN c fi 2 /tc 2 is the density of states at the Fermi surface fi. The explicit form of the 
function F(y) is 



POO 

F(y) = y 2 dx 



x 



y/x 2 - y 2 (e x + l) 2 



where y = (3 A = A/ksT and x = Ep/ksT. The gap A depends on the temperature. In limiting 
situations, the temperature dependence of the superfluid gap A for 1 Sq pairing is given by 



A(T) = A(0) - (27rfc B TA(0)) 1/2 exp 



A(0) ' 
k B T 



3.06 k B T c (1 - T/T c 



,1/2 



I (T < T c ) 
(T c - T < T c 



(19) 



with A(0) = 1.76kBT c . In lieu of this, F(y) ~ (1 — ^) as T — > T c so that it vanishes quadratically 
with A. For T <C T c , F(y) vanishes exponentially as e _2A//T . Accurate fits of F(y) valid over the entire 
range of y may be found in Ref. ||. 

It is instructive to compare the result Eq. ( |TTD with that of Flowers et al. |T| for nonrelativistic 
neutrons: 



(GfCvX 



N ,^_L,; (0W (_|_ 



F(v) 



1 + H(£1) 2 (^' J 

42 VCy/ V C 



(20) 



where z/(0) = rn*pF n /n 2 is the density of states, m* is the Landau effective mass 0, and -UF n = PF n /fn* 
is the velocity at the Fermi surface. Since the PBF process occurs in the vicinity of the Fermi surface, 
the formal similarity of the results between the nonrelativistic and extreme relativistic cases is not 
surprising. The neutral current couplings (see Table 1), the density of states at the Fermi surface, 
and the results of phase space integrations are, however, somewhat different between the two cases. 
For both relativistic quarks and nonrelativistic neutrons, the dominating contribution arises from the 
vector channel although the suppression of the axial channel is more severe in the nonrelativistic case. 

The main approximation made in the emissivity calculation for massless quarks is that A <C fi. 
Modifications due to a finite quark mass m enter in the dispersion relation e p = y/p 2 + m 2 — \jp 2 F + m 2 , 

with a corresponding density of states z/(0) = N fN c pp^J p 2 F + m 2 /7i 2 , and from corrections to the quark 
matrix element through the spinor normalization as Y> s u s (p)u s (p) = j$ + m. In what follows, we set the 
spin along the direction of the momentum of the massive quark. The quark field operator in the spin 
basis may be written as 



V(x) 



E 
p 



e ip - x u p {xp,+bp,+ + Xp,~h 



-} + e lp -%{XvA, 



XR-4-} 



(21) 



3 Defined at the Fermi surface, m* accounts for momentum-dependent interactions in a potential model 
description. In a field-theoretical description, m* accounts for scalar interactions that off-set the bare mass. 
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This acts on a quasiparticle state \p + ,p'_) and restores it to the condensate. 
The vector channel squared matrix element is 



M_(p) = U-(p)u-{p) 



) 2 • c 2 

p 1 + m 



{\M q \ 2 )^ = 2(u p v p , + v p u pl f ■ c 2 v Tr[M^pf)YM.(p)f 



1 - 75^ (p) 



(22) 



with n(q) = ^qj- Setting m = reproduces the chirality projector that was used for the cal- 
culation in the massless case. Taking the nonrelativistic limit m ^> |p|, we find that only /i = = v 
contributes to the matrix element. This is consistent with the nonrelativistic matrix elements consid- 
ered in the neutron superfluid case Q. 

The squared matrix element in the axial channel is 



2(u p Vp> — VpUpi) 2 



c 2 A Tr[M_{p/)^ l5 M„{p)n 5 } ■ 



(23) 



To perform the trace, notice that there are 0(m 2 ),0(m), and O(m ) terms that may be grouped 
separately. The 0(m) terms in the vector and axial channels vanish, because they involve the trace of 
an odd number of 7-matrices or are anti-symmetric in the very indices in which the overall neutrino 
integral is symmetric. 

In the vector channel, the 0(m 2 ) term is 



m 



—{4^ + 4n^n' u + An u n'^ - A(n.n')g^} 



(24) 



while the O(m ) term is 
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(25) 



The axial channel gives similar contributions as above, but with the following changes: (i) the first and 
last terms in Eq. fl24|) are both of opposite signs, and (ii) the last term in Eq. (p5|) is of opposite sign. 

Simple expressions for the emissivities that take into account quark masses may be obtained by 
retaining leading terms in some limiting situations. 



(a) For m/pF 1, the projector in Eq. (|22|) reads |(1 + 75) + y. In this case, the total emissivity is 



e R (m) 



e R (Q) 



m 

4p| / 



1 fc A 



7m 2 



(26) 



where e R (0) is the result in Eq. (|17D . In the context of the long-term cooling of neutron stars with 
quark cores, the utility of this result lies in the extent to which it applies for massive strange quarks. 



(b) For p F /m <^ 1, the projector in Eq. (|22| ) becomes m(l + 70) (1 + 757 • p) 
In this case, 



IpI(7-P-7o7s)- 



e v (m 



e A {m) 



ey(0) 



m 



31 



P 2 F 



cv 



Tn 2 + p 2 F 84 m 2 + p 2 F 

v F \ 2 11 fm 2 + p 2 F /2\ 
~) 42 \m 2 +p 2 F J 



and 



P 2 f 



m 2 + pp 



(27) 



where z/(0) and F(y) in ey (0) involve the dispersion relation for massive particles. The changes from the 
extreme nonrelativistic result in Eq. (|20f) are easily understood. The factor m 2 +p 2 F in the denominator 
of the two terms inside the square bracket comes from energy factors in the Lorentz invariant phase 
space, while the numerator receives contributions from 0(m 2 ) and O(m ) terms in the trace. The 
heirarchy A < < m should be maintained for these expressions to be valid. We can express the 
sum of the emissivities from the vector and axial channels in Eq. (p7|) as 



e NR (m) = e y (0) 




(28) 



In this form, we recognize the nonrelativistic result of Flowers et al. in Eq. (p0|) , and in the axial 
channel an additional term involving (m*/m) 2 , where m* is the Landau effective mass (see footnote 
2) defined through p F = m*vp/c = (vr 2 ^) 1 / 3 . The ellipsis denotes terms that are smaller by further 
powers of (pp/m) 2 compared to the preceding term in Eq. (|28f) . The presence of the (m*/m) 2 term 
and its importance in the case of superfluid protons for which cy is vanishingly small, was first noted 



in Ref. [12] and was obtained as a relativistic correction in PBF studies of a proton superfluid [12 



For further discussions on additional corrections, see Ref. ||13|| . 

The results in Eqs. (|17|) and (|2~6D are the prinicipal new results of this work. For simplicity and for 
numerical comparison with emissivities from other processes in quark and baryonic matter, we will use 
the result in Eq. fllTT) in what follows to represent the vv emissivity from the PBF process in quark 
matter. In cgs units, Eq. ( p]) may be expressed as 



ef = 1.4 x 10 20 N v TlF a Q ( —) erg cm" 3 s" 1 , (29) 

\n J 



where N u is the number of neutrino flavors, T 9 is the temperature in 10 °K, F is the result in Eq. 

a q = c v 1 + | (j^) (see Table 1), and n B = p Fq /ir 2 is the baryon density (n = 0.16 fm~ 3 is 

the empirical nuclear matter equilibrium density). This may be compared with the emissivities in 
normal quark matter from the bremsstrahlung (q + q — > q + q + v + v) and modified Urea (e.g., 
d + d—*d + u + e~ + v) processes [§] 



brems „ 3 x 1Q 19 ^8 _B ^ ^-3 g -l ^ (30) 



s brems _ 3 x 1Q 19 ] y rpS ( ^V'' 

q \n J 
ef Urca ~ a 2 c cos 2 6 C e b g rems , (31) 

where a c is the strong interaction coupling and 6c is the Cabbibo angle. Note that both of these 
processes are suppressed by superfluidity by a factor of ~ exp(— 2A//cbT). Thus, as T falls below T c , 
the PBF process dominates the cooling in quark matter. 

In baryonic matter, the vv emissivity from the PBF process is 

e 



? = 1.65 x 10 21 ^T 9 7 ( ^) 1/3 [^Ftai Y^ 3 — \ erg cm" 3 s" 1 , (32) 



n J \-r m 



where = c 2 ri and Yi = rii/nB is the baryon concentration of species i. The values of the neutral current 
hyperon couplings in Table 1, which were calculated in [14j by considering both the SU(3) singlet and 
octet contributions to the baryon current operators, suggest that in matter containing hyperons, the 
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PBF process can occur with hyperons of all species. The values quoted in Table 1 of |13|], where the 



A and S° couplings were found to vanish and the S~ coupling was very small, differ significantly from 
those in Table 1 probably due to the neglect of the singlet contributions in Ref. [|1^]. It is intriguing 
that at the same baryon density in baryonic and quark matter, e B v /e u q v ~ (constant) • (rig /no) -1 / 3 hes 
in the range 2 — 5 for the densities of interest in the long-term cooling of solar mass compact stars. 

The bremsstrahlung (n + n — >n + n + u + u) and modified Urea process (n + n — > n + p + e~ + v) 
emissivities are 

i. , x in (^) (^) 1/3 I?,,,, (Tg cm-'' s 



\m n J \n 

e MUrca ^ 1Q 22 ( V * erg cm -3 g -l (33) 

V n J 

where R nn is a density dependent reduction factor, that lies in the range 0.3-0.15 for densities in the 
range 0.5 — 3 no, recently computed in Ref. || utilizing on-shell nn amplitudes (similar reductions are 
expected for the modified Urea process also). 

So far, we have restricted ourselves to the case of isotropic 1 S pairing, which provides a larger 
gap compared to anisotropic p-wave pairing ( 3 P 2 states). The fact that the smallest gap is the most 
relevant for long-term cooling highlights the importance of p-wave pairing though it requires higher 
energies (or densities) than s-wave pairing. For example, for neutrons in the lowest energy mj = 
quantum number of the 3 P2 state, the anisotropy of the gap yields the relation || 

A(T) = 1.98k B T c (l — T/T c ) 1/2 (T c - T < T c ) . (34) 

Just as for nucleons, the quark field operator in Eq. (^T|) must be modified to matrix-valued creation 
and annihilation operators due to the triplet nature of pairing. It can be obtained by a unitary 
transformation on the ground state [jnj with the coherence factors in Eq. ([|) acquiring additional 
factors from the angular dependence. The p-wave emissivity exhibits a T-dependence that is similar 
to that from s-wave pairing. 

The effects of strong and electromagnetic correlations on weak interaction processes involving 



baryons are known to be large [TfJ. By considering only the electromagnetic correlations between 
protons and electrons, Leinson has recently demonstrated that vv emissivity from the PBF process 
from gapped protons is substantially larger than that obtained at the mean field level JL7||. A similar 
enhancement may be expected in quark matter through a collective response of u, d, and s quarks in 
a charge neutral environment even with miniscule concentrations of electrons. 



Calculations to include effects germane to "QCD condensed matter" |TD| would have to address 
modifications in the weak interaction hamiltonian that will likely lead to couplings and quasiparticle 
excitations that are modified from those considered here. We defer calculations of PBF processes in 
such phases and detailed cooling simulations of stars with quark matter to subsequent studies. 

Research support from the DOE Grant No. FG02-88ER-40388 is gratefully acknowledged. We 
thank Dima Yakovlev for useful communications. 
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TABLE 1: The standard model neutral current vector and axial- vector couplings of neutrinos to 
quarks and baryons. The hyperon couplings are taken from Ref. [I3| and include both the octet and 
singlet contributions, t stands for e— , or r— type neutrinos. Numerical values are quoted using 

D=0.756, F=0.477, and sin 2 0^=0.23. 
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-\{F + § ) = -0.3645 


V£ + £~ - 


->■ + £~ 


-| + 2 sin 2 W 


= -1.04 


±(L> - 3F) = -0.34 


VI + T.+ - 


-» i/< + £+ 


\ — 2 sin 2 # w 


= 0.04 


+ = 0.6165 


vt + S° - 


-> !/< + S° 


l 

2 




\(D-F) = 0.14 


vi + H" - 


-> ^ + H~ 


-| + 2sin 2 


= -1.04 


±(.D - 3F) = -0.34 


^ + s°- 


-i/< + H° 


l 

2 




-±(£> + F) = -0.6165 
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